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ABSTRACT 

Recently, it has been shown that soft-state black hole X-ray binaries and active galactic 
nuclei populate a plane in the space defined by the black hole mass, accretion rate 
and characteristic frequency. We show that this plane can be extended to hard-state 
objects if one allows a constant offset for the frequencies in the soft and the hard state. 
During a state transition the frequencies rapidly move from one scaling to the other 
depending on an additional parameter, possibly the disk-fraction. The relationship 
between frequency, mass and accretion rate can be further extended by including 
weakly accreting neutron stars. We explore if the lower kHz QPOs of neutron stars 
and the dwarf nova oscillations of white dwarfs can be included as well and discuss 
the physical implications of the found correlation. 
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1 INTRODUCTION 

The general idea that X-ray binaries (XRBs) and active 
galactic nuclei (AGN) have similar central engines, if scaled 
with black hole mass, is increasingly supported by empirical 
correlations connecting both classes. First, the fundamen- 
tal plane of accreting black holes connects XRBs and AGN 
through a plane in the black hole mass, radio and X-r ay lu- 
minosi tv space ([Merloni et al. . 2003; F alcke et al.ll20o3 ). Re- 
cently, iMcHardv et al.l (120061 ) have reported that a second 
plane connecting stellar and supermassive black holes exists 
in the space defined by the accretion rate, the black hole 
mass and a characteristic timescale of the X-ray variability. 
This suggests that all black holes can be unified by taking 
the accretion rate as well as the black hole mass into account. 
Here, we explore the latter correlation in further detail and 
discuss the inclusion of strongly sub-Eddington black holes 
(BHs) and neutron stars (NSs). 

The power spectral density (PSD) of NS and BH 
XRBs can be well described by a nu mber of Lorentzians 
with variable coher ence factor Q (e.g., iPsaltis et ail 1 19991 : 
iBelloni etT al. 2002b). Each Lorentzian is described by their 
characteristic frequency (the frequency of the maximum of 
the Lore ntzian in the f r equenc y times power plot). As pre- 
sented in IBelloni et all (|2002bh . it is possible to fit NSs and 
BHs with 

• a zero-centered low frequency Lorentzian Lit fitting the 
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low-frequency end of the band-limited noise with character- 
istic frequency v^. Usually, this frequency is denoted as vt 
in the literature. However, to avoid confusion with the break 
frequency studied in AGN, we add the prefix "1" for low. 

• two Lorentzians fitting the high-frequency end of the 
band-limited noise with frequencies vi and v u (the lower 
and upper high-frequency Lorentzian Li and L u ). These two 
Lorentzians take over the role of the upper and lower kHz 
QPO if they are present. 

• one or two Lorentzians fitting the region around the 
frequency of the low-frequency quasi-periodic oscillation 
(QPO). The narrow core of the QPO (Llf) has the charac- 
teristic frequency vlf, while the broader "hump" Lorentzian 
Lh has a characteristic frequency Vh- 

For a sketch showing the different Lorentizans see Fig. Q] 

XRBs are observed in several states: the hard, 
the soft and the intermediate state (I MS). For defin i- 
tions, subcla s ses a n d examples see e . g.. | Nowakl (|l995h : 
Belloni et all (120051): THoman fc Bellonil (|2005l ), but see also 



McClintock fc Remillardl (120061 ) for slightly different defini- 



tions. 

The measured shapes of AGN PSDs are similar to those 
of sof t-state XRBs (see e.g. lMcHardvlll988l : lMarkowitz etafl 
I2003D . The timescales for AGN are several orders of mag- 
nitudes longer, roughly in agreement with the expectation 
that all length-scales increase linearly with black hole mass. 
The uncertainties of the measured AGN PSDs are usu- 
ally too large to fit a detailed Lorentzian model (but see 
McHardy et al. in prep.), due to the long timescales. Thus, 
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Figure 1. Sketch of the different Lorentzians found in BH and 
NS XRBs. The x-axis denotes the frequency and on the j/-axis 
we show frequency times power. Besides the Lorentzians found in 
hard-state objects we also show a broken power law model found 
for the PSDs of soft-state BHs. We have increased the normali- 
sation of the soft-state model so that it lies above the model for 
the hard state. 



one uses broken power law models. The shape of the PSD 
of AGN and XRBs at high Fourier frequencie s is well de- 
scribed by a power la w with index ~ — 2 (e.g.. iGreen et all 
ll99llCui et alj|l997h . For Seyfert Galaxies and those soft- 
state XRBs similar to Cyg X-l the PSD at low Fourier fre- 
quencies turns to a power law with index ~ — 1 (see Fig. [T 



lEdelson fc Nandral 19991 ; lUttlev et al ]|2002l ; ICui etalTl99 
The Fourier frequency of the turnover will be denoted as the 
high-break frequency v^. Using th e bolometric lum i nosity 
as a tracer for the accretion rate, iMcHardv et al.l (|2006h 
found that this high-break frequency is related to the ac- 
cretion rate and black hole mass as VhbM oc M/Mvdd for 
soft-state objects. 

It has recently been found bv lMigliari et all (|2005l ) that 
for the BH GX 339-4 the frequency of the lower high- 
frequency Lorentzian v\ and for the NS 4U 1728-34 the 
frequency of the "hump" -Lorentzian Vh is correlated with 
the radio luminosity in their hard state as i/^j, oc L R 7 , 
where v^h denotes either vi or v^. As the radio luminosity 
is a good tracer of the accretion rate (|Kording et al.ll2006l . 
Lr oc M 1 ' 4 ), the correlation can be written as vi : h oc M. It 
is therefore suggestive that the "variability plane" found for 
high state XRBs and AGN can be extended to hard-state 
BHs or even NSs. 

As shown in Fig. [T] the superposition of the different 
Lorentzians (Lit, Lh, Li and L u ) creates a flat plateau in the 
vP v plots, and the PSD can therefore be roughly described 
with power law with index ~ — 1 for frequencies between i/u, 
and vi. This power law steepens to a power law with index 
~ — 2 at the peak of the high-frequency Lorentzians (Li and 
L u ) similar to the behaviour found in soft-state XRBs and 
AGN. It is therefore likely, that the frequency of either the 
upper or lower high-frequency Lorentzian (v u or v%) corre- 
spond to the high break frequency found in soft-state XRBs 
and AGN. The lower high-frequency Lorentzian L; usually 
dominates the high-frequency end of the band limited noise 



(e.g.. IPottschmidt et al.ll2Q03 | ). Additionally, v i and v u are 
not linearly coupled (e.g., Belloni et al.1 l2002bl ). While the 
frequency of the lower high-frequency Lorentizan vi seems 
to be proportional to M , v u is therefore not linearly pro- 
portional to M for individual sources. This suggests, that 
the lower high-frequency Lorentzian L; with the character- 
istic frequency v\ corresponds to the high break frequency 
of AGN. 

We note that for NS one observes "para llel tracks" 
for kHz QPOs on freque ncy flux diagrams (e.g.. iFbrd et al.l 
2000; Ivan der Klisl koOl ,. ie., while flux and frequencies are 
strongly correlated for short periods of time, the correlation 
is nearly absent for longer timescales. If the X-ray luminos- 
ity (flux) is a good tracer of the accretion rate, there cannot 
be a simple one to one dependence of the frequency vi of the 
lower high-frequency Lorentzian on accretion rate. However, 
there are at least three possible tracers of the accretion rate 
(frequencies, X-ray luminosity and radio luminosity) and the 
frequencies or the radio luminosity might be better indica- 
tors _Lh^ii_the_2^ayljirminosjtY ; 

iPsaltis et all \ 19991 ); iBelloni et all (l2002bh present a 
correlation of the frequencies of the lower kHz QPO with 
those of the low frequency QPO vlf- This correlation can be 
extended to include noise features of BHXRBs by using the 
lower high-frequency Lorentzian L; as the lower kHz QPO. 
Also cataclysmic variables (CVs) c an be included by u sing 
dwarf nova oscillations (DNOs, e.g. JWarner et al"1l2003T ) and 
the QPO frequency. Here, we explore if the varability plane 
found in sort state AGN and XRBs can also be extended to 
hard-state BH XRBs and NSs or even to white dwarfs. 



2 THE SAMPLE 

To explore the dependence of the variability properties on 
the accretion rate and mass of the compact object, we con- 
struct a sample of hard and soft-state BHs and NSs with 
estimated accretion rates. 



2.1 Accretion rates 

Soft state BHs and all NSs are generally assumed to be effi- 
ciently accreting. Thus, we can use the bolometric luminos- 
ity directly as a measure of the accretion rate: 

Lboi 



M = 



0.1c 2 : 



(1) 



where we assumed an accretion efficiency of r\ = 0.1. 

For hard-state BHs the accretion rate is not linearly 
related to the X-ray as the acc retion flow is likely to be 
ineffic ient (e.g.. lEsin et al.l Il997l ). However. IKording et al.l 
(2006) provide estimates of the accretion rates from either 
radio or X-ray luminosities. The accretion measure based on 
the radio luminosity is 



M = 3 x 10 



12/17 



(2) 

Vl0 30 erg e- 1 J s y ' 

In the notation given in IKording et al. I (|2006h we have set 
/ = 1 and rj = 0.1 to allow for a direct comparison between 
hard-state BHs and sources with a measured accretion rate 
obtained from the bolometric luminosity. Using the funda- 
mental plane of accreting black holes (and the radio/X-ray 
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correlation for XRBs) the radio luminosity accretion mea- 
sure can be translated to an accretion rate estimate based 
on the 2-10 keV X-ray luminosity: 

0.43 



M « 3.4 x 10 1 



/ L 2 



( M 



(3) 



Vl0 36 ergs-V \Max339 , 

We note, that this accretion measure does not only depend 
on the X-ray luminosity but also on the mass of the black 
hole. While the masses of BH XRBs are typically around 
~ IOM0 this mass-term provides an additional uncertainty 
compared to the accretion rate measure based on the radio 
luminosity. Thus, we will estimate the accretion rate from 
the radio luminosity if quasi-simultaneous measurements of 
the radio flux and the timing features are available. Other- 
wise, we have to use the X-ray flux to obtain an accretion 
rate estimate. 

The uncertainty of these accretion rate measures is 
hard to access as there are only a few data-points available 
to normalise the accretion measure. The sample standard- 
deviation of the 4 BH data-points is ~ 0.3 dex, but has 
little significance. If one includes also the NS points (overall 
14 points) we find a scatter of ~ 0.2 dex. Another possibility 
to measure the uncertainty of the accretion measure is via 
the fundamental plane of accreting BHs. If both accretion 
measures (based on the radio and X-ray emission) presented 
here are exact, there would be no scatter in radio/X-ray 
correlation found for XRBs and the fundamental plane for 
XRBs and AGN. The scatter around the fundamental plane 
can therefore be used as rough estimator of the uncertainty 
of the accretion rate measure. The cleanest sample for the 
fundamental plane including only hard-state XRBs and low 
lumin osity AGN has a scatter of 0.15 dex (|Kording et al.l 
2006). This suggests that the uncertainties of the accretion 
rate measure based on radio luminosity is around ~ 0.2 dex. 
The accretion rate measure based on the X-ray luminosity 
will have a similar intrinsic uncertainty, but has additional 
uncertainties due to uncertainties of the BH mass estimates. 

The accretion rate measure based on the radio lumi- 
nosity is also applicable to island state NS. However, as the 
radio luminosity accretion rate measure has been normalised 
using accretion rates obtained from bolometric luminosities 
its absolute accuracy is unlikely to exceed the measure based 
on bolometric luminosities for efficiently accretion objects. 
In summary we will use the following accretion rate mea- 
sures in order of preference: 

(i) Accretion rates from bolometric luminosity for soft- 
state BHs and NSs 

(ii) Accretion rates estimated from the radio luminosity 
for hard-state BHs 

(iii) Accretion rates estimated from the X-ray luminosity 
for hard-state BHs 

Throughout this paper we will measure the accretion 
rate in g/s, all frequencies in Hz and masses of compact 
objects in solar masses. 



2.2 Black Holes 

Our hard-state BH sample is based on measurements in the 
public RXTE archive as well as already published data. We 
select all sources which have well measurable frequencies v\ 
as well as estimates for the black hole mass and distance. 



However, in order to obtain an estimate of Vi one needs a 
hard-state power spectrum, since in the hard-intermediate 
state characteristic frequencies are higher and the broad 
lower high-frequency Lorentzian (Lj) cannot usually be de- 
tected (with the exception of GRS 1915+105). Moreover, 
the source count rate needs to be sufficiently high to allow 
a detection. This rules out the final parts of the outbursts 
and limits our sample to bright early hard-state observa- 
tions. Unfortunately, only few sources have been observed 
sufficiently early in their outburst to fulfill all the require- 
ments: 



• GX 339-4: IMigliari et~al1 (|2005l ) present values for v x 
and the radio fluxes of GX 339-4 in its hard state. These 
radio fluxes have been used to estimate the accretion rate. 
The d istan ce to GX 339-4 is still unc ertain, IShahbaz et al.l 
l|200lf) and lJonker fc Nelemansl (|2004l ) give a lo wer limit of 6 
kpc, but the distance may be as high as 15 kpc (|Hvnes et al.l 
l2004h . We therefore ado pt a distance of 8 kpc. The mass 
function is 5.8 ±0.5 Mq (jHvnes et aJj|20 03j), which is there- 
fore a lower limit for the mass of the black hole. If one as- 
sumes zero mass for the companion and a mean inclina- 
tion angle we obtain a mass of M « 12Mq. Additionally, 
this higher mass fits the fundamental plane as well as tim- 
ing correlations better than smaller masses. Thus, we use 
M=12M e . 

• XTE J1118+480: \Bellom et ail (|2002bh measured v x for 
XTE J1118+480 on the 4th and 15th May 2000. The VLA 
has observed the source during that time. We have analysed 
archival 8.5 GHz data and found a flux of 5.7 mjy on the 
27th April, 6.4 mjy on the 13th May, and 6.8 mjy on the 
31st May. We interpolate linearly between those dates and 
estimate the 8.5 GHz flux on the 4th of May to be 5.9 mjy 
and 6.4 mjy on the 15th May and estimate the accretion 
rate fro m these radio fluxes . We assume a mass of 6.8 ± 
0.3 Mr7, ilRitter fc Kolbll2003t ) and a distance of 1.71 ± 0.05 
kpc (|Chatv et al.ll20031 ). 

• XTE J1550-564: For the outburst of XTE J1550- 
564 in 2002 we use vi values and X-ray fluxes from 
Bcllo ni. Colombo. Homan. Campana fc van der Klisl 



(2002a). This outburst never left the hard state. Ad- 
ditionally, we include the outbursts of 1998 and 2000. 
For these outbursts, we reanalysed the data of the RXTE 
archive to obtain vi. The 2-10 keV X-ray flux has been 
estimated from the 2-9 keV PCA counts assuming that 
the X-ray spectrum can be described by a power-law with 
photon index of T = 1.5. As the observed X-ray band is 
similar to the energy range of the estimated flux, the spec- 
tral uncertainties do not strongly effec t the estimated flux. 
We a ssume a distance of 5.3 ± 2.3 kpc Jjonker fc Nelemansl 
2004h and a mass of 10.6 ± 1 M Q (|Orosz et al.ll2002l ). 

• GRO 31655-40: For GRO J1655-40 we reduced public 
RXTE data on MJD 53429.7 to measure vi and obtained 
the 2-10 keV flux from the spectral model given on the 
webpag^B mai ntained by J. Hom an. We assume a mass of 
6.3 ± 0.5 Mm ([Greene et alj|200lb and a distance of 3.2 kpc 
|Hiellming fc Rupenlll995l '). but see lFoellmi et al.1 l|2006l) for 
significantly lower distance estimates (<1.7 kpc). We there- 
fore assume 3.2 ± 1.5 kpc. 



1 http : //tahti .mit . edu/opensource/1655/ 
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• GS 1354-6 U: For GS 1354-644 we assume a mass of 
7.34 ± 0.5 M dCasares et al.1 12004) The distance to the 
sourc e is uncertain, [ Citamoto et alT(|l990l ) suggests 10 kpc 
while ICasares et al.l ( 2004 ) give a lower limit of 27 kpc. We 
will assume a distance of 10 kpc. The 2-10 keV flux and the 
measurement of vi have been obtained from public RXTE 
data. 

• XTE J 1650-500 RXTE started to observe the source 
during the outburst in 2001 just as the source starts it tran - 
sition from the hard to the soft state l|Homan et al.| [2003). 
We use their first observation which is still in the hard state 
to obtain a value for vi and the 2- 10 keV flux. The dist ance 
to J1650-500 is be tween 2-6 kpc dTomsick et all 120031 '). we 
assume 4 ± 2 kpc. lOrosz et al.l l|2004h suggest that the BH 
mass of J1650-500 is between 4 and 7.3 Mq, we use a mass 
of 5.5 ± 2 M_0. 

• Cyg X-l: To compare these hard-state measurements 
with a soft-state object, we also show Cyg X-l in its 
soft state. The high break frequencies are taken from 
lAxelsson et al.l l|2006r ). To obtai n accretion r a tes w e use the 
bolometric luminosities given in IWilms et al](|2006h , as soft- 
state objects are likely to be efficient accretors. The avail- 
ab le data has been binne d in luminosity bins as described 
in iMcHardv et alj f2006!). We only include the source for 
comparison, as the model describing the PSDs for the soft 
state is different to the Lorentzians we use here. As black 
hole ma ss we use IOMq, an d the distance is assumed to be 
2.1 kpc (|Massev et al.lll995l ). 

• GRS 1915+105: This source is always accreting 
at a large fraction o f the Eddington accretion rate 
l|Fender fe Belloni I200I ). It is therefore likely to be effi- 
ciently accreting, so we can use th e bolometric luminos - 
ity to estimate th e accretion rate dKording et all [2006). 
We o btain vi from iBelloni et al.l (|2002bh and X-rav fluxes 



from iTrudolvubovl (|200ll) . We assume a ma ss of 15 M ( 



and assume the distance to be 11 kpc ( Fender et al.l 
1999l:lDhawan et al.ll2000l: IZdziarski et alJlibod . but see also 
Chapuis fc Corb el 2004; Kais er et al.ll2004l for lower values) 



2.3 Neutron Stars 

For all NSs we assume a mass of 1.4 Mq. 

• low luminosity X-ray burster. We include IE 1724- 
304, GS 1826 - 24, SL X 1735-269 with v x measurements of 
IBelloni et ail (l2002bl ). The 2-10 keV X-ray flux has been 
obtained from archival RXTE PCA data. To obtain the 
bolometric luminosity and therefore the accretion rate, we 
use a bolometric correc tion of 2.5 found for atoll NSs 
l|Migliari fc Fenderl 120061 ) also for these X-ra y burster. We 
assum e distances of 6.6 kpc f or IE 1724-304 dBarbuv et al.l 
1 19981 ). 5 kpc for GS 18 26-24 jThompson et alj 120051 ) and 8 
kpc for SLX 1735-269 jMolkov et al.ll2005T ). 

• IGR J00291+5934: The accreting millisecond X-ray 
pulsar IGR J002 9 1+5934 went into outburst in 2004 
<|Shaw et al.ll2005l ). iLinares et ail (|2006l ) have measured its 
frequency 1/1, and we obtain 2-10 keV flux measurements 
from archival RXTE data. As V507 Cas is also in the field 
of view of the PCA, we assume that the quiescent flux can 
be attributed to the V507 Cas and subtract this from the 
measured flux of J0 0291+5934 during out burst. We assume 
a distance of 4 kpc ( Ga lloway et al.ll2005l ). 



• atoll sources: For the atoll s ources we include 4U 1608 - 
522 with vi measure ments from van Straaten etail (|2003h 
(distance of 3.4 kpc; Ijonker fc Nelemansl l2004h . The 2-10 
keV flux has been estimated from the 3-9 keV PCA count 
rate assuming a power law with F = 1.6. To obtain a bolo- 
metric luminosity we use a bolometric correction of 2.5. For 
4U 1812-12 we use public RXTE data to obtai n vi measure- 
ments and use the bolometric X-ray fluxes from lBarret et ail 
i|2003t ). It should be no ted that our measur ed values for ui 
are double of those of Barret et al. (|2003l ). We assume a 



distance of 4 kpc (|Cocchi et al. 200C ). 

Z-sources: We include the Z-sources GX 340 



Jonker et al.l 19981. distance of 9.5 kpc) and GX 5-1 
Wiinands et al. I ll998l . 7.4 kpc). The conversion factor from 



the 2-16 keV PCA count rate to bolometric luminosity 
of our Z-s ources has be e n obt ained from the Z-source 
GX 17+2 |Pi Salvo etail <|2000l ). The 0.1-200 keV flux 
given in the paper has been divided by the 2-16 keV PCA 
counts. For the Z-sour ces we use the lower k Hz QPO as the 

frequency vi following IBelloni et al. (l2002bh. 

• 4U 1728-34: For 4U 1728-34 iMigliari et ail i|2005l ) mea- 
sure the dependence of vu on the radio luminosity and find 
that it scales with vu cx L°ji ad oc M. We include these 



pints for comparison . We assume a distance of 4.( 
Galloway et~aill2003l ). 



kpc 



2.4 AGN 

For AGN, we use the sample presented by 
lUttlev fc McHardvl (120051 ). They give the high break 
frequencies, as well as the black hole masses and the 
bolometric luminosities. The majority of sources are Seyfert 
galaxies, and thus they are likely high- or very-high- 
state objects. We can therefore translate the bolometric 
luminosities directly to accretion rates. 



3 RESULTS 

For our A GN and the two stella r soft-state black holes 
we found in IMcHardv et all (|2006l ) that the measured high 
break frequencies u^t lie on a plane in the i^b> M and M 
space: 



log v h6 = £ acc log M + £ m log M + b v , 



(4) 



where £ acc and £ m denote the correlation indexes for the 
accretion rate and the black hole mass; b v denotes the 
constant offset. We found that face = 0.98 ± 0.15 and 
f m = —2.1 ± 0.15. Both parameters are within the uncer- 
tainties of integer value. We will therefore adopt the integer 
solution (f acc = l,f;m = —2) and will not present new fits 
here. 

The edge-on projection of the plane using our sample 
is shown in Fig. [3] The two lines indicating Vi^M oc MM~ 2 
{^i,ht denotes vi for hard-state objects and Vhb for soft-state 
objects). The upper line is the fit to soft-state stellar BHs 
and supermassive BHs using integer valued parameters for 
the variability plane. For the adopted parameters and our 
units the constant offset b v for the soft-state objects is b v — 
— 14.5 ± 0.1. This offs et is different from the one given in 
IMcHardv et~ail l|2006h . as we use different units and fixed 
the parameters of the plane to integer values. 
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Figure 2. Left side: Our sample of stellar BHs. For Cyg X-l we plot v^, while we show vi for the other objects. The lines indicate 
Vi oc MM~ 2 . While the upper line is a fit to soft-state objects (including AGN, see iMcHardv et al. (2006)), the lower is a fit to the 
hard-state XRBs only. Right side: Same as left side, but here we plot stellar and supermassive objects. In both panels soft-state and IMS 
objects are plotted with open symbols. 




f5.5 16.0 16.5 17.0 17.5 18.0 18.5 19.0 19.5 ?5.5 16.0 16.5 17.0 17.5 18.0 18.5 19.0 19.5 

Log M [g/s] Log M [g/s] 

Figure 3. Left panel: NSs with measured vi in comparison with our stellar BH sample. The NS sample seems to follow the correlation, 
albeit with larger scatter. The upper line shows the correlation normalized to soft-state objects while the lower line is normalized for 
hard-state objects. Right panel: Inclusion of kHz QPOs of Z sources. While the kHz QPOs of Z sources are near the correlation this is 
not true for kHz QPOs of atoll sources. See text for a detailed discussion. 



3.1 Stellar black holes 

On the left side of Fig. [2] we show the projection of the 
variability plane "zoomed-in" on our sample of stellar BHs. 
For hard-state BHs we have argued that the frequency vi 
of the lower high-frequency Lorentzian corresponds to the 
high-frequency break Vhb used for soft state. In Fig. [2] we 
also plot our hard-state BH sample. Similar to the scaling 
found in soft-state XRBs and AGN, we find that our hard- 
state BHs follow a scaling vi oc M . However, the constant 
offset is different for hard and soft-state objects, and we also 
show a plane normalized to the hard-state objects only. The 
constant offset between the soft and the hard state scaling 
is —0.9 dex. 

The uncertainties of the measured values shown in 
Fig. [2] are dominated by systematic errors (e.g., uncertain- 



ties in the accretion rate measure) as well as uncertainties 
of the black hole mass and distance of the black holes. Most 
uncertainties of all measurements of one source are cou- 
pled. Showing these coupled uncertainties as errorbars for 
the data-points is therefore misleading. The error budget 
consists of: 

• Uncertainties of the primary parameters Vi and the ra- 
dio or X-ray flux: These uncertainties are different for each 
data-point. A typical value is ^ 0.04 dex. 

• Uncertainties in the BH mass and distance measure- 
ments: These uncertainties are shared for all data-points of 
a given source. Typical values are between 0.03 and 0.3 dex. 
This component dominates the total error budget, especially 
as they enter the variability plane quadratically and our two 
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main sources (GX 339-4 and XTE J1550-564) have large un- 
certainties in the mass estimate (0.2 and 0.15 dex). 

• Uncertainties of the accretion rate measure: This un- 
certainty has two components: Firstly, an uncertainty of the 
normalisation of the accretion rate measure, this would di- 
rectly change the normalisation of all data-points. Secondly, 
uncertainties due to source peculiarities and the individual 
measurements. The exact value of these uncertainties is hard 
to access. In Sect. 12.11 we have argued that their combined 
effect is roughly ~ 0.2 dex. 

The correlation of ui oc M is found in both hard-state 
BH that have several measurement (GX 339-4 and XTE 
J1550-564). The majority of the uncertainties just men- 
tioned do not play a role if one only considers only a single 
source as they would only change the normalisation constant 
b v . Thus, we can safely assume that there is a correlation 
between vi and M for individual sources. Whether the con- 
stant offset b v is significantly different for hard and soft-state 
objects need to be verified. 

To access the value and uncertainty of the constant off- 
set b v for our sample of sources we use two approaches: First 
we can compute the direct mean and sample variance of the 
sources around the correlation using the same weight for 
all data-points. We find b v = —15.39 ± 0.04 with a sample 
variance of 0.18. However, as many of the errors of the data- 
points are coupled we might underestimate the uncertainty 
of the mean. 

To avoid the coupling of the uncertainties for different 
data-points of a given source we first treat every source indi- 
vidually and combine the different sources in a second step 
taking the estimated errors into account. To calculate the 
error on the mean b v we used the following steps: 

(i) Calculate b v for each source individually and measure 
the intrinsic scatter around the linear correlation through 
the sample variance (pint)- 

(ii) Estimate the uncertainty of the measured b v for all 
sources. We set 



c& = Cint + (2ctd) 2 + (2<tm) 2 + 



(5) 



where go is the uncertainty of the measurement of the dis- 
tance, om the mass uncertainty and the systematic un- 
certainty of the Af-measure. 

(iii) Calculate that weighted mean and error of mean with 
the estimated uncertainties. Each source is weighted by their 
uncertainties and the number of measurements of the source. 

Using this method, we find b v — —15.38 ± 0.08. The mean 
value of b v is in agreement with the value found using the 
sample variance. Finally, we have to consider that the over- 
all normalisation of the accretion rate measure has an un- 
certainty of ~ 0.2 dex. This uncertainty is shared by all 
data-points and all errors in the estimate of the normali- 
sation effects the mean value of b v directly. If we include 
this uncertainty we obtain an overall error of 0.22 dex. As 
the difference between the normalisation found for the hard 
and the soft state is 0.9 dex, we conclude that the difference 
between the hard and the soft-state scaling is significant. 

Using this meathod we can also calculate the expected 
mean variance of the correlation given the assumed uncer- 
tainties. We find a — 0.38, which is larger than the sample 
variance of 0.18. The discrepance is mainly due to the cou- 



pled uncertainties, but we may also have overestimated the 
uncertainties of the accretion estimator or those of the dis- 
tance and mass measurements. 

Only for GX 339-4 and XTR J1118+480 we have used 
the radio luminosity to obtain accretion rates while for the 
other hard-state objects we used the 2-10 keV flux. We have 
verified that Fig. does not change significantly if we use 
the accretion rate measure based on the X-ray flux for all 
sources. This is not surprising, as the accretion rate measure 
used for the X-ray flux is partly based on the radio/X-ray 
correlation of GX 339-4. Thus, our findings are not affected 
by our choice of different the accretion rate estimators. 

The outbursts of XTE J1550-564 in 1998 and 2000 seem 
to follow a different track than the linear dependence seen 
in the hard-state outburst in 2002 and the other hard-state 
objects. The first two outbursts rise more steeply than the 
linear hard-state scaling. The outbursts seem to start at a 
frequency vi characteristic for a hard state and lead to val- 
ues typically found for soft states. In these outbursts the 
photon index changes strongly during the observations. In 
1998 the firs t observation in our sample has a photon index 
of T = 1.53 (|Sobczak et alJl200Ch which changes continuesly 
to the last observation which has F = 1.98, which is atypical 
for the hard state were one would expect the photon index to 
stay hard until the state transition. In the 2000 outburst the 
source is also softening from F = 1.46 to 1.7 while the power 
law cutoff moves from 33 keV to 19 keV with a significant 
black-body component visible in the last observation used 
(a typical cut-off for the hard state is ~ 80 keV). Thus, it is 
likely that the source was already in an intermediate state 
in these outbursts. This suggests that there is one scaling 
relation for the classical hard state and one with a higher 
normalisation for the soft state, but there is obviously an- 
other parameter that governs the transition between those 
scaling relations. During the transition from the hard to the 
soft state, the source increases its frequency fast and moves 
from the hard scaling to the soft one. 

For XRBs the transition from the soft to the hard scal- 
ing is probably best traced by the spectral index F. A strong 
correlation between timing properties and the photon in- 
dex has been found in objec t s nea r the intermediate states , 
see e.g.. iPottschmidt et all <|2003h : IVignarca et~aH (|2003h : 
iKalemci et all (|2005l ). A simple dependence on the photon 
index is however not directly applicable to AGN. The tim- 
ing properties are also c orrelated to the hardness ratio see 
e.g.. iBelloni et a l. (2005). The hardness ratio can be gener- 
alised to the non-thermal fraction f — T — Lp K , where 

J L Di3k+ L PL ' 

L/Disk is the luminosity in the multi-color black body com- 
ponent while Lpl describes the luminosity of the power law 
component. This non-thermal fraction has similar proper- 
ties as the hardness rat io but is also applicable to AGN 
(see iKording et all 12006'). Sources solidly in the soft state 
(/ = 0) follow a track given by log v l>hh = log (MM -2 ) -14.7, 
while those in the hard state (/ = 1) follow log^.^j, = 



log(MM~ 



14.7 — 0.9. In the transition between both 



states, the frequencies depend on the photon index and 
therefore on /. This suggests that the dependence of v on 
M, M and / can be approximated by: 



log n„ 



log M - 2 log M - 14.7 - 0.9 8(f) 



(6) 



Where 9 is a monotonic function with 8(0) — and 8(1) = 
1. Unfortunately, with the readily available data of state 
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transitions we are not able to fully constrain it. We note 
that this dependence may have to be further modified due to 
hysteresis effects similar to those seen in hardness-intensity 
diagrams. 

3.2 Neutron Stars 

In the left panel of Fig. [3] we show our sample of NSs with 
visible broad lower high-frequency Lorentzians (L;) together 
with the stellar BHs. The low luminosity X-ray burster and 
the accreting millisecond pulsar lie near the expected scaling 
for hard-state BHs - but with increased scatter. Also our 
atoll sources are in a similar frequency range. However, while 
we could observe the linear dependency of vi on accretion 
rate in some single sources for BHs (ie., XTE J1550-564 and 
GX 339-4) this is harder for the NSs, since the accretion rate 
changes in our sample are small. The only object changing 
its accretion rate significantly is IGR J00291+5934. For that 
source, the slope seems to be slightly shallower than what is 
found for the BHs. However, as already mentioned, a second 
source (V507 Cas) is in the field of view of the PCA, so it 
is hard to measure correct fluxes for the pulsar. It may be 
that we are underestimating the lowest detected fluxes. As 
a sample the NS sources still scatter around the correlation 
foun d for hard-st a te BH s, albeit with a larger scatter. 

IPsaltis et al.l (|l999h present a tight correlation between 
the frequencies of the low frequency QPO (vlf) and those 
of the l ower kHz QPO in NS s and BHs. In the unified pic- 
ture of iBelloni et~ai1 l|2002bl ') the lower kHz QPO found in 
some NSs is identified with the broad lower high-frequency 
Lorentzian L\ for BHs and weakly accreting NSs. As the 
correlation between vi and MM~ 2 holds for some sources 
on the vlf/vi correlation, it may be possible that one can 
extend our accretio n rate vi correlation to all sources on the 
IPsaltis et ail (l999) correlation. 

We have just shown, that atoll sources with a measured 
frequency vi roughly follow the variability plane. At higher 
accretion rates atoll sources show a lower kHz QPO at sig- 
nificantly higher frequencies than the observed values of ui 
at low accretion rates (when the lower kHz QPO is observed 
the broad lower high-frequency Lorentzian Li is not seen). 
However, the increase in the accretion rate is not sufficient 
to ensure that the lower kHz QPO is in agreement with the 
correlation. This can be also shown by looking at Vh for 1728- 
34. Vh is the frequency corresponding to the broad "hump" 
Lorentzian (Lh) of the low frequency QPO. Thus and the 
frequency of the low frequency QPO are well correlated and 
have roughly similar values. As we show on the right side of 
Fig. [3] Vh is (at least for 1728) near the correlation. Thus, 
the lower kHz QPO will be above the correlation by a factor 
~ 15. 

Z-sources have significantly higher accretion rates than 
atoll sources and often show kHz QPOs. The frequencies of 
the lower kHz QPOs of Z-sources are shown in Fig. [3] The 
accretion rates are so large that the frequencies of the lower 
kHz QPO of Z-sources are near the hard-state correlation. 
However, Z-sources are very strongly accreting, e.g., they 
belong to a bright IMS. These sources should therefore fol- 
low the high state slope similar to the BH GRS 1915+105. 
The lower kHz QPO frequencies are so high, that they have 
the same order of magnitude as the Keplerian frequency at 
the surface of the neutron star (~ 1700Hz). Thus, it may 
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Figure 4. Edge-on projection of the variability plane including 
AGN, stellar BHs and NSs with directly measured 

be that the frequencies start to saturate and are therefore 
slightly below the high state correlation. See also the discus- 
sion on WD and the corresponding Fig. HJ] where this effect 
may be even more pronounced. 

The slope of the Z-sources alone seems to be steeper 
than the correlation found for BHs. It may be, that while 
the accretion rate does not change very much the sources 
makes a state transition from the analog of the hard IMS 
to the soft IMS (horizontal branch to normal branch). As 
we have seen for BHs, during state transitions the source 
moves rapidly form the hard-state scaling to the soft-state 
scaling. This rapid transition explain the steeper slope of Z- 
sources. Nevertheless, as the kHz QPOs of atoll sources do 
not fit on the correlation we can not be sure if it is not pure 
coincidence that the Z-sources seem to follow the correlation. 
For the correlation we therefore have to rely on sources with 
measured vi. 

We have mentioned in the introduction that the "par- 
allel tracks" found in flux frequency diagrams may suggest 
that the frequency of the kHz QPO may be a better tracer 
of the instantaneous accr etion rate than the X-ray lumi- 
nosity (|van der Klisll200ll L We are currently unable to use 
any frequency as a direct tracer of the accretion rate as the 
conversion factors are yet unknown. As we have shown it 
is likely that the frequency (ui) does not only depend on 
the accretion rate but also on the state of the source. Thus, 
at least vi cannot be used as a measure of the accretion 
rate for all times if one does not have detailed information 
on the state of the source. If the frequency of the lower 
high-frequency Lorentzian does not only depend on the ac- 
cretion rate but also on an additional parameter describing 
the accretion state, one could explain the parallel tracks. On 
short timescales the main parameter is either the accretion 
rate or the accretion rate and the accretion state parameter 
are strongly coupled creating the well correlated tracks. On 
longer timescales the relationship of the two parameter can 
change so that one observes parallel tracks. 

The final edge-on projection of the variability plane in- 
cluding BH XRBs, NS XRBs with directly measured vi and 
AGN is shown in Fig. [4] As a sample the NSs extend the 
correlation toward lower accretion rates. 
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Figure 5. Edge-on projection of the variability plane compared 
to our sample of WDs. The uncertainty on the accretion rates in 
a factor 2-3, so in logspacc 0.3-0.5 dex. 



Table X. Parameters of our WD sample. 





M [g/s] 


"UNO Hz 


log uM 2 


(log uM 2 ) pred 


IX Vel 


5 x 10 17 


3.7 xlO" 2 


-1.62 


3 


UX UMa 


4 x 10 17 


3.4 xl0~ 2 


-1.56 


2.9 


OY Car 


3 x 10 16 


5.5 xlO" 2 


-1.59 


1.8 


VW Hyi 


> 6 x 10 14 


~lx 10" 2 


-2.1 


0.1 



Masses: IX Vel: 0.8 M©, UX UMa: 0.9 Mp, OY Car: 0.69 M ( 
and V W Hyi: 0.84 Afo. Refs: IX V el: iBeuermann fc Thomas 
(199C), UX UMa: iFroning et al.l ll2~003h; ISuleimanov et al 
[ 2004): iBaptista et al.1 <1998fl; iKnigge et all (Il998h. OY Car 
iMarsh fc Home! dl998l); iPratt et al.1 ll99E )iTwood et al.l l|l989|) 



Woudtl il2006l) : IPandel et alj 



2003) 



VW Hyi: IWarner 
ISion et alj ll 19971) 



3.3 Comments on DNOs for accreting white 
dwarfs 

The aforementioned correlation of the fr equencies of the LF 
QPOs with those of the lower kHz QPO (IPsaltis et al.lll999l') 
has a lso been extended to include CVs (e.g.. IWarner et aU 
2003). Here, the dwarf-nova oscillation (DNO) frequency 
corresponds to the frequency vi of the lower high-frequency 
Lorentzian. In table Q] we tabulate the parameters of four 
well known systems. The accretion rates of CVs have higher 
uncertainties than those of NSs or BHs, as the majority of 
the bolometric luminosity is not emitted in the well access- 
able X-ray band. We assume that the accretion rates are 
only correct up to a factor of 2-3. Especially the accretion 
rate for VW Hyi is uncertain, as we use the quiescent value 
as a lower limit to the accretion rate, while the DNO fre- 
quencies are measured during the end of the decline from an 
outburst. Thus, the accretion rate will be approximately the 
quiescent rate. Nevertheless, observed values of log vM 2 are 
mostly around -1.6, while we expect values between 0.1 and 
3 (assuming that they follow the soft-state scaling). Thus, 
in the current form, the correlation is not valid for CVs (see 
Fig. [5]). We note that, for CVs, the DNO frequency rises 



with accretion rate (|Woudt fc Warner! liobl ). similar to the 
expected behaviour of the correlation. 

However, up to now we have not considered that a white 
dwarf is significantly larger than a neutron star or a black 
hole. It is not evident how to correct the correlation for the 
larger radius. For BHs the radius is linearly related to the 
BH mass. Thus, the correlation found for BHs suggests that 
one or both mass terms of the correlation correspond to a 
radius. If one of the mass terms corresponds to a radius, 
the correlation reads vi cx j^, where R is the radius of 
the central object. This correlation can be read as a corre- 
lation between vi and the power liberated in the accretion 
disk around the central object divided by M 2 . The power 
liberated in a standard accretion disk around a central ob- 
ject with radius R is R/(6GM/c 2 ) times smaller than the 
power liberated in a disk around a non-rotating BH, the one 
radius term corrects for this fact. A white dwarf of average 
mass has a radius roughly ~ 500 times larger than that of 
a NS. This would move the expected position of OY Car 
and VW Hyi near our correlations (see Fig. [5} • Their posi- 
tions are in agreement with the correlation given the large 
uncertainties on the accretion rates. IX Vel and UX UMa, 
however, are still far from the correlation. Those two sources 
are nova-likes, with high accretion rates. If the correlation 
would hold for all those sources, their frequencies would ex- 
ceed the breakup frequency of the white dwarf (~ 5) Hz. If 
the DNO frequency (or v{) is somehow related to the Ke- 
plerian rotation this is impossible. Thus, it is likely that 
the sources follow the correlation until the frequencies are 
comparable to the fastest possible frequency (breakup fre- 
quency for WDs and the frequency of the innermost stable 
orbit for BHs). For these high accretion rates the frequen- 
cies can not increase any further and stay constant around 
the fastest possible frequency. Schematic frequency tracks 
for WDs, BHs and NSs are shown in Fig. [U] The shown 
lines should only illustrate possible frequency tracks. As the 
break up frequencies are reached at high luminosities the 
shown tracks are normalised to the high state correlation. 
The shown function is a smooth broken power law which is 
linear with M at low accretion rates and constant at high 
accretion rates (y cx M(l + (M/Mn) 2 )~ 0,5 ). If one assumes 
that both mass terms of the correlation correspond to radii, 
none of the measured WDs fall on the expected correlation. 
Thus, we consider it to be more likely that only one mass 
term corresponds to a radius. In summary, while dwarf no- 
vae systems seem to lie near the correlation if one takes the 
size of the central object into account, nova-like objects do 
not seem to follow the correlation as their frequencies would 
exceed the break up frequency of the WD. 



4 DISCUSSION 

In the previous section we found that the variability plane 
found in soft-state accreting black holes (Mc Hardy et al.l 
2006) can be extended to hard-state stellar black holes (see 
Fig-01- The frequencies depend on the accretion rate (M) 
and the black hole mass (M) as: 
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Figure 6. Projection of the variability plane onto the v\ -M^n 
plane. Besides the scaling relation expected for hard-state objects 
we also show the breakup frequency of WDs (~ 5 Hz) and NSs 
(~1700 Hz) and the frequency at the innermost stable orbit of 
BHs (~ 240 Hz). Sources can only follow the correlation until the 
frequencies are comparable to the breakup frequency, there the 
correlation need to level off to the constant breakup frequency. 

for soft-state objects; for hard-state objects the proportion- 
ality factor is ~ 0.3. It is even possible to extend the correla- 
tion to NS with directly measured frequencies vi (see Fig. [4}. 
Here, we discuss the physical implications of these findings. 

This dependence on the accretion rate as well as on the 
black hole mass can be simplified by as suming that the ac- 
creting black hole is scale invariant fe.g.. lMerloni et al. 2003; 
iFalcke et ai]|2004r ). If the frequency arises through Keplerian 
motion or any other process scaling with the size of the cen- 
tral object, the left side of the equation is invariant of black 
hole mass. The right side denotes a simple linear dependence 
on the accretion rate in Eddington units. 

While a linear dependency on the accretion rate is seen 
if a source is in the canonical hard or soft state this is not 
the case for sources in intermediate states or near the state 
transition. For example, we have seen that XTE J1550-564 
follows the linear scaling in its hard outburst in 2002, while 
the frequencies rise faster than linear for the other outbursts. 
Thus, besides accretion rate there must be at least a sec- 
ond parameter governing the observed frequencies (vi,v hb ). 
Similar to the hysteresis found for state changes, (see e.g, 
the hardness intensity diagrams of XRB outbursts,) there 
se ems to be a hystere sis effect for the frequencies (see Fig. 7 
in iBelloni et afll2005l ) . 

In case the high-frequency break originates from Keple- 
rian motion or a frequency that is related to the Keplerian 
frequency, we can use eq. (O to obtain the dependence of 
the active radius r a , which corresponds to the break in the 
PSD, on the accretion rate. This yields: 

/ s -2/3 

= 625 ( ) , (8) 
Rg \0MM Edd J 

where Rg is the gravitational radius Rg = GM/c?. The 
correlation with the accretion rate therefore implies that 
the radius associated with the break moves inwards with 
M -2 / 3 . Interestingly, such a behaviour may be observed in 



GRS 1915+105. IBelloni et all (|l997n reports that the disk 
radius associated with the region emitting the measured X- 
ray spectrum changes from 70 km to roughly 20 km while the 
RXTE/PCA count rate increases by a factor of ~ 5. If the 
bolometric correction does not change significantly between 
the different observations, this suggests that the emitting 
radius moves inwards as r -0,75 compared to our prediction 
ofr" - 67 . 

It is tempting to associate the active region creating 
the high-frequency break with the transition from the geo- 
metrically thin standard disk to the optically thin inefficient 
accretion flow. However, most of our objects are soft-state 
objects. Within the soft state it is usually assumed that the 
optically thick, geometrically thin accretion disk reaches up 
to the last stable orbit. The luminosities of a given source in 
the soft state can change by a factor 15 (2-30% Eddington), 
so the inner disk radius would have to change by a factor 6. 
This is not in agreement with the redii deduced from X-ray 
spectral fitting nor with observed scaling of the X-ray lumi- 
nosit ies with disc temperature as T 4 re.g jGierliriski fc Pond 
2004). Furthermore, the observed frequencies are far too low 
to originate from a Keplerian motion at the inner edge of 
the accretion disk (innermost stable orbit). We further note 
that for a normal standard disk, viscous timescales are too 
long (~ 2 s) to explain the high frequencies and seem to have 
a very low dependency on accretion rate (a _4 ^ 5 M _3/ ' 10 ). It 
may well be that the break frequency is indeed associated 
with a Keplerian motion. However, as we have argued, the 
corresponding radius does not seem to be directly associated 
with the inner edge of the standard disk. Maybe this radius 
only reaches the innermost radii for highly accreting objects 
like Z-sources with their lower kHz QPO. 

We have seen that NSs with measured v\ seem to fol- 
low the suggested correlation. At higher accretion rate the 
Z-sources have lower frequencies than expected. This dis- 
crepancy may arise from the slightly larger size of the NS 
with its boundary layer compared to a BH with similar mass. 
While increasing the accretion rate, the emitting region can- 
not move inwards as expected as it is stopped by the bound- 
ary layer. This may also be the reason for the relatively 
low DNO frequencies of strongly accreting white dwarfs (see 
Fig. El). 



5 CONCLUSION 

We have shown that the variability plane found in AGN 
and soft-state XRBs (y-^bM oc M/Mem) can be extended 
to hard-state BH XRBs and weakly accreting NSs. For 
hard-state objects the frequency of the lower high-frequency 
Lorentzian (yi) corresponds to the high-frequency break in 
AGN and soft-state XRBs. For white dwarfs, the situation 
seems to be more complicated as we need to correct of the 
WD radius. While some dwarf novae seem to lie near the 
correlation (within the large uncertainties), strongly accret- 
ing nova-like systems do not seem to follow the correlation 
as the frequencies would otherwise exceed the breakup fre- 
quency of the WD. In case that the high-break frequency 
in the PSD is related to a Keplerian orbit around the black 
hole, we find that this radius in gravitational radii depends 
on the accretion rate as (M/MEdd)' _2,/3 ^ However, the orbit 
corresponding to the high break is unlikely directly related 
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to the radius where the standard optically thick disk turns 
into an optically thin inefficient accretion flow. 
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